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Abstract: We quantify perturbative and non-perturbative QCD effects in the ex- 
clusive J/V-photoproduction cross section, and in the shrinkage of the differential 
cross section with respect to momentum transfer, t. We predict that in the high 
energy THERA region there will always be a significant contribution to this process 
that rises quickly with energy. This implies that the taming of the rise of the cross 
section with energy, due to both the expansion of spatially-small fluctuations in the 
photon and to higher twist effects, is rather gradual. 
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1. Introduction 

One of the major challenges facing particle physics is to understand the interplay 
of perturbative and non-perturbative QCD effects in the high-energy production of 
heavy quark bound states and in the interaction of such states with hadrons. Under- 
standing the details of this interplay is interesting in its own right because it reveals 
the practical boundaries of applicability of perturbation theory. It is also necessary 
to achieve an unambiguous interpretation of many related phenomena, including the 
suppression of the yield of J/ip mesons produced in nucleus-nucleus collisions. The 
investigation of hard exclusive processes, such as the exclusive photoproduction of the 
J/i[) meson in photon-proton collisions considered here, gives a unique opportunity 
to quantify this physics. 

For several years, high energy, hard exclusive and semi-inclusive processes, which 
include photo- and electroproduction of heavy vector mesons and deep inelastic pro- 
duction of light vector mesons and real photons (DVCS), have been modelled reason- 
ably successfully using the exchange of two gluons in a colour singlet configuration 

|, & |, § |, g §. 



For vector meson production initiated by longitudinally polarised photons || Q 
and for DVCS (Kj such calculations have now been placed on a firmer theoretical 
footing by the proof of QCD factorization theorems which show that perturbative 
two gluon exchange is the dominant process in the asymptotic limit (Q 2 — > oo). 

To make a process "hard" it is necessary to provide a large momentum scale 
(either a heavy quark mass or Q 2 , or both) which squeezes the hadronic fluctuation 
of the photon so that small perturbative qq configurations are responsible for the 
dominant contribution. For the diffractive photoproduction of J/ip this hard scale is 
thought to be provided by the charm mass. However, since the charm quark is rather 
light (m c p» 1.5 GeV), and transverse polarisations of the quasi-real photon domi- 
nate, there are likely to be significant contributions from non-perturbative regions in 
which the hadronic fluctuation of the photon has a large transverse size (these be- 
come progressively less important as the photon virtuality increases). This interplay 
between soft and hard contributions has been seen in a phenomenological way by 
the success of the two-Pomeron fit of Donnachie and Landshoff |ll] in which the soft 
Pomeron term, associated with non-perturbative effects, appears make a significant 
contribution to J/ip photoproduction at HERA energies. 

At very high energies, or small x Bj , the time taken for a given fluctuation in the 
photon to interact with the proton target is considerably smaller than its formation 
time and the time required to produce the hadronic final state. This implies univer- 
sality for the interaction cross section over a wide range of inclusive and exclusive 



processes. In |~2| we introduced a model for this universal cross section, a, for all 
transverse size fluctuations*. We produced a satisfactory description of the inclusive 
cross section data, indicating that our model is reasonable. In this paper, having 
made some suitable minor adjustments, we apply the model to J/ip photoproduc- 
tion. Our aim is twofold, to provide a good description of this process (and hence 
to elucidate the role of non-perturbative physics in J/ip photoproduction) and to 
further constrain the model in order to be able to make better predictions for other 
exclusive processes in the future, using the same framework. 

Our model for the interaction cross-section is based on the well-known leading- 
log perturbative QCD result for the interaction of a small transverse-size qq dipole 
which proceeds via two-gluon exchange ]T3[: 



° P QCD(b ,x) = —b a s (Q ) xg(x', Q ) , (1.1) 

where scales x' and Q 2 , which depend on transverse size b, are described in Subsection 



37T| . This form is applicable for transverse sizes b < 6qo ~ 0.4 fm. For larger trans- 



*Modelling is currently unavoidable in QCD due to the necessity to take non-perturbative con- 
finement effects into account. We defer a more detailed modelling of the strong interaction in the 
region of small a s , which seems to be required for the very small- x regime in which Structure 
Functions may achieve the unitarity limit, for future studies. 



verse sizes we introduced an ansatz based on the known behaviour of soft hadronic 
interactions (we matched on to the measured np cross section at b = b^ = 0.65 fm) 
and introduced a suitable interpolation for 6q < b < b^ and an extrapolation for 
b > b n . In this paper we slightly simplify the model for a by connecting the points 
6qo and b^ with a straight line (instead of extrapolating a fit to the known 6-shape 
just below 6q as i n 13 )• At ver y n i§ n energies, we tame the steep increase due 
to the rapid rise of the small-x gluon density by imposing a unitarity restriction, 
a(bi,W 2 ) = (r(b 7T ,W 2 )/2, and connect the point 61 < 6qo to b^ using a straight line 
(see figs.(^,|l^) and Subsection |3]4] for more details). 



We use this model, based on eq.Ql.lj) with CTEQ4L [pi gluon input density, 
evolved using skewed evolution, to investigate the interplay of perturbative and non- 
perturbative effects in J /if) photoproduction^. At the higher photon-proton energies 
of the HERA range (W ~ 250 GeV) this involves probing the gluon distribution at 
very small x ~ 10 -4 , i.e. outside of the range in which it has been tested directly. In 
this region, it is mainly constrained indirectly by its effect on the observed scaling 



violations of the inclusive structure function F 2} predicted by DGLAP |TJ] evolution. 
It should be remembered that the predictions we make would change if the input 
gluon density is changed. Indeed, in our final plot for the cross section we also 
show the results using the latest MRST leading order partons [TB|] : ' : . Taken together, 
the predictions using both CTEQ4L and MRST produce a spread which spans the 
currently available data. Other models for the interaction cross section, a, may be 
found in the literature (see e.g. [|18 , |i~9l and references therein). 



The amplitude for exclusive processes involves a convolution of a with light-cone 
wavefunctions for the initial-state photon (known from QED [^0| in the qq case) 



and for the final-state diflractively-produced object. For the light-cone wavefunc- 
tion of the J /if) we use a hybrid wavefunction introduced in ||. This is derived 



by solving the Schrodinger equation using the Buchmueller-Tye potential model pi 
(constituent quark mass m c = 1.48 GeV), boosting the resulting Schrodinger wave- 
function to a fast moving frame. The wave function at small b < 0.3 fm is fixed 
by imposing QCD behaviour (<fiy(z,b = 0) oc z(l — z)) and normalized using the 
known leptonic decay rate of J/tf). In this way we account for the QCD radiative 
corrections to the qq component of the charmonium wavefunction for small trans- 
verse sizes b < 0.3 fm. This effectively takes into account the strong modification of 
non-relativistic charmonium-model predictions for the leptonic width due to QCD 
radiative corrections (for a recent review see ||22|| ). The exclusive formation of the 



heavy bound state, the details of which are embodied in the light-cone wavefunction, 



^The value of <J p QCD{b 2 ,x) should not change significantly within the next-to-leading order 
approximation, as compared to the leading order, because it is determined from fitting the same 
value of i*2- However, both the relative size of <j p qcd(P 2 ', x) and xg(x, Q 2 ) and the numerical value 
of xg(x, Q 2 ) itself may change rather significantly. 

■^This leading order input distribution is more recent than, but similar to, that found in [p~7J. 



suppresses the contribution of higher order Fock states (such as \ccg>) in the virtual 
photon, making this a particularly good process in which to test the universality of 
the dipole picture. 

For small enough x, the rapid rise of a due to the small-x rise of the gluon den- 
sity should be tamed to avoid violating the unitarity restriction for the interaction 
of spatially-small partonic fluctuations of the photon. Two types of effects act to 
tame the increase of high energy processes. As the energy increases, spatially-small 
partonic fluctuations in the photon typically expand rapidly in transverse size due 
to the increased phase-space for radiation. This expansion in size, known as Gribov 
diffusion, is a consequence of the randomness of radiation. For large size configu- 
rations perturbative QCD is inapplicable so the amplitudes no longer have this fast 
increase with energy. We account for Gribov diffusion to some extent by implement- 
ing a behaviour typical of soft hadronic cross sections (i.e. with only mild increase 
with energy) at large transverse sizes. Even if a is completely independent of b 2 
at sufficiently small-x the Structure Function of a hadronic target should retain a 
residual ln(l/x) increase with energy as a result of the infinite normalization of the 
wavefunction of a virtual photon (the infinite renormalization of electric charge is 
due to hadronic vacuum polarisation). Another effect is taming, or shadowing, of the 
rapid rise with energy due to higher twist effects. So, for small enough x, unitarity 



corrections were introduced in |L2| to tame the rapid rise of a due to the small- x 
rise of the gluon distribution. These affect larger transverse sizes in the perturbative 
domain first. If perturbative QCD effects are tamed by the unitarity of S-matrix 
only, this mechanism would lead to an increase of the typical impact parameters, p, 
involved in the scattering of qq dipole from a nucleon as p 2 oc In 2 1/x and leads to a 
related increase of Structure Functions as F 2 oc (In 1/x) 3 (and a' oc In 2 1/x). 

In this paper, we illustrate that the photoproduction of the relatively light J/ijj 
can act as a precursor for this taming since it is sensitive to relatively large trans- 
verse sizes, not just to the wavefunction at the origin j|, ||. Thus, we consider the 
behaviour of the energy dependence of the cross section of J / -^-photoproduction (es- 
pecially that of the slope of the momentum transfer, t, which is parameterised by 
a'jp) to be crucial in establishing the existence of a new regime in which the standard 
DGLAP approximation is violated. Remember that the slowing down of the increase 
of the Structure Function alone does not necessarily imply violation of DGLAP, 
because it may be due to Gribov diffusion, which is a leading twist phenomenon. 
We note in passing that an important role of large transverse distance effects (soft 
QCD) reveals itself in the energy dependence of the slope for the photoproduction 
of J/ip mesons which, within our model, should be intermediate between the soft 
regime, ct^p(soft) ~ 0.25 GeV -2 , and the perturbative regime, a'p <C ct^p(soft) (see 
Subsection [T3] and Section 6). 

We predict a reduction in the steepness of the energy dependence of amplitudes of 
hard exclusive processes which will begin to take effect at the higher HERA energies 



and in the THERA region (250 < W > 1000 GeV). The precise details are of course 
specific features of our model, which incorporates a simple taming ansatz in the small 
x region. More generally, taming corrections are eventually expected to reduce the 
rise with energy of all hard small-x cross sections at given impact parameter §. 

This paper extends the work of [pi, || in two original directions. Firstly, it ex- 
plicitly includes the non-perturbative component coming from large transverse sizes 
and provides a reasonable unified description of the process from the low energies 
measured at fixed target experiments [^3j as well as the HERA data [^4[ [25[ ^6f| . 
Within the framework of our analysis the relative contribution of non-perturbative 
effects may be quantized, albeit in an inevitably model-dependent way. Secondly, 
we make predictions for energies beyond the HERA range W > 300 GeV that may 
eventually be tested at a higher energy ep collider such as the proposed Tesla-HERA, 
or THERA, project (see e.g. f27P ). The logic of |12| dictates that taming corrections 



are required in this high energy region and we qualify and quantify the expected 
effect of these corrections on J/ip photoproduction. 

The paper is organized as follows. Section 2 contains the basic formula for the 
cross section. Section 3 investigates various issues surrounding the implementation 
of this formula including setting scales in the gluon distribution, running quark mass 
effects, t— dependence, skewedness and the calculation of the real part of the ampli- 
tude. In Section 4 we illustrate the effect of changing the rescaling parameter A, in 
the model for the dipole cross section. Section 4 also contains our predictions for the 
THERA energy range. Following a discussion and evaluation of a'jp in Section 5, we 
include a general discussion in Section 6 and conclude in Section 7. 

2. Basic Formulae in the photoproduction limit 

From eq.(50) of [||] in the limit Q 2 —>■ 0, neglecting the real part of the amplitude, 
the differential cross section for the photoproduction of J/ip reads: 



da , 12n 3 T v M v 



|t=o 



dt a e . m (4m, 



c 2 ) 4 



7T 2 



[^) 2 C(Q 2 = 0), (2.1) 



where My,Yy are the mass and leptonic decay width of the vector meson and 
m c = 1.5 GeV is the charm quark mass. Using the notation of || (in particular 
eqs. (25,26,51)) we have for the overall dimensionless suppression factor 

C(Q 2 = 0) = ( ^^ ] T(p)R(0) , (2.2) 



^However, it is unclear whether the unitarity of ^-matrix for the interaction of small size fluctu- 
ations should tame the increase with energy of Structure Functions, because of the related increase 
of the important impact parameters with energy. This question will be investigated in a separate 
publication. 
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where m Cjr (b) is the running charm quark mass and (jE->4>v are the light-cone wave- 
functions for the transversely-polarised photon and vector meson, respectively. The 
latter depend on transverse size b, and on z, the momentum fraction of the photon 
energy carried by the quark. In eq. ( |2.1|) the gluon density and a s have been extracted 
at the average point, <b>, of the integration over b in the amplitude, using the re- 
lationship Q 2 a = A/ (b 2 ), with A ~ 10. In |J we attempted to go further than this 
average approximation by sampling these functions at the Q 2 = A/6 2 , underneath 
the integral in transverse size b. The factor m\ r implicitly depends on b, so rightfully 
also belongs underneath the integral in b, sampled at an appropriate scale. 

Reinstating a underneath the integral according to this procedure leads to a 
modified version of eq.(|2.1|): 
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where the penultimate line makes use of eq.(40) of |J. The imaginary part of the 
amplitude, Sm A, and its normalization, N, are given by 
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where the light-cone wavefunction for the photon, 



K (b m Ct r), is purely trans- 



verse for photoproduction. In this equation we have chosen to separate out a piece 
of the 6-integral, I z (b), which only depends on the light-cone wavefunctions of the 
vector meson and photon and is independent of energy. 

For small dipole sizes b < 0.3 fm QCD behaviour is imposed on <pj/f, so it is 
appropriate to run the charm quark mass underneath the b-integral too, in both the 



overall m 2 r factor and in the argument of the Bessel function, using the appropriate 
renormalization group equation for masses (see Subsection p.2|) . 
Finally, assuming the usual exponential fall-off in t we have 

J(7 P - m P) = ^ + .^" A? (29) 

where the real part of the amplitude has been reinstated via /3 = 3?e A/Ssm A 
(see Subsection ft.5Q . The HI collaboration recently reported [Q a value of B — 
4.73 ± 0.25 Zq 39 GeV~ 2 . Recently ZEUS reported [^| an improved measurement 
of B in J/ip photoproduction and found it to depend on energy. We will incorporate 
this shrinkage using a simple form (see Subsection 



3. Improvements to the basic formula: uncertainties in the 
cross section 

In this section we expand and explain various features and uncertainties in the basic 
formulae of Section 2. Many crucial issues involve the 6-integral in eq. (|2.7| ). Firstly 
we explain the choice of scales, x', Q 2 used in the dipole cross section and make a first 
comparison to the available data, using very basic assumptions. Then we consider 
the changes induced by considering a more careful treatment of running mass effects, 
shrinkage, non-zero real part and skewedness. The latter requires that we replace the 
ordinary gluon with the skewed gluon in eqs.( |l.l|]2~77| ). We follow the usual choice of 
conventional input distributions evolved using skewed evolution. 

Lastly, in the next Section, we decrease the parameter A which relates transverse 
sizes to four-momentum scales, via Q 2 = A/6 2 , from 10 to 4. This has the effect of 
increasing the non-perturbative contribution to the cross section and leads to a much 
better description of the current data. We also give predictions for THERA for both 
values of A. 

3.1 Scale setting in the gluon density 

In IJT2J we examined inclusive structure functions, Fl and F 2 , using our model for 



a. We employed the following prescription for the b-dependence of the scales x' and 
Q 2 : 

x' = <*»(! + 0-75^^) , (3.1) 

Q 2 = ^, (3.2) 

with, 
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By examining the standard perturbative QCD formula for Fl for light quarks 
we found < x >~ 1.75 x (we designed the ansatz to reproduce this for large Q 2 
when b = (b)). We also found that a value of A = 10 reproduced the perturbative 
QCD results fairly well. However the results turned out to be rather insensitive to 
the precise value of A. In this paper we explicitly examine the sensitivity of J/ip 
photoproduction on A and find that it is much more sensitive to it than Fl and F 2 
(see Section 4). This is interesting because the value of A determines the dividing 
line between perturbative and non-perturbative physics. 

In the photoproduction limit, Q 2 — ► 0, for charm eqs.( [3.1| , |3.3|) lead to 
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for the b— dependent momentum fraction of the incoming gluon. This will be used 
for the sampling of a underneath the b- integral in eq. (|2.7|) . 
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Figure 1: A comparison of the J/ip photoproduction cross section, using conventional 
CTEQ4L parton distribution function (PDF) for the gluon and j3 = 0, with data f2|, |||, 
25, 26]. The solid curve has fixed mass and slope B = 4.0 GeV~ 2 , the long and short 
dashed curves include running charm quark mass and TF-dependent slope, respectively. 
The dotted curve includes both effects. 



The solid line in fig.([l]) shows the resulting J/ip photoproduction cross section 
using eq.( |2.9| ), with CTEQ4L gluon density, (3 — 0, a fixed slope parameter of B = 4.0 
GeV -2 and fixed quark mass of 1.5 GeV. The available data |23|, [24], ^, |26| is also 
shown. This curve undershoots the data at low energies and overshoots at high 



energies, so it appears to rise too steeply with energy to provide a good description 
of all available data. We now consider several features which improve the shape of 
the energy dependence. 

3.2 Running quark mass 

In || hard QCD corrections were introduced in the vector meson wavefunction: for 
b < b (b = 0.3 fm was chosen for J/ip) QCD-behaviour was imposed (0y oc z{l — z)) 
and normalized to the measured leptonic decay width. For consistency it was then 
necessary to replace the constituent quark mass with the running quark mass for 
small b. This was implemented in the following way: (cf. eq.(38) of ||) 



ia 



In this paper, we demand that the mass satisfies the renormalization group 
formula for quark masses (see e.g p8|): 



mcAQ ] = mc U(Sj) J ' ( } 

where Qf ~ 2.0 GeV corresponds to the matching scale, b Q = 0.3 fm, at which 
hard corrections are applied. We set the number of flavours rif = 4 both in the 
exponent and in the one-loop beta-function of a s (choosing n/ = 3 instead would 
make very little difference). This effect suppresses the small b region relative to 
the case in which the mass is taken as fixed, imposing the desired QCD behaviour 
m c — > for Q 2 — > oo. Fig.(||) shows m 2 cr {b)/m 2 c versus b. We use m c = 1.5 GeV, 
which is approximately the constituent quark mass of Buchmueller-Tye potential 
model Pl| , which we use to construct </y. Fig.(^|) illustrates the size of the effect 



at the amplitude level. Running the quark mass influences both the argument of 
the Kq Bessel function and the overall multiplicative factor in eq. Q2.8p . It turns out 
that there is only significant suppression for b < 0.2 fm (rather than b < bo = 0.3 
fm). A comparison of the solid and long-dashed lines of fig.(|T|) illustrates the rather 
small decrease that this change induces in the cross section as a result of the reduced 
contribution from small b. 

3.3 Energy dependent slope 

We shall evaluate this effect in detail in the section 5. Here we give a brief outline of 
the significance of the energy dependence of the slope of the ^-dependence. The recent 
presentation of ZEUS data, from the 1996-97 running period, at DIS2000 |25], p6| 
indicated evidence for shrinkage. Following a fit to da/dt of the form exp(Bt) in 
each of the seven bins in W, a value of a'p = 0.098 ± 0.035 ± 0.05 GeV -2 is obtained 
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Figure 2: The ratio of running charm quark mass squared to fixed mass squared as a 
function of transverse size using eq.( |3.5| ) and incorporating four light flavours, rif = 4. 



from examining the energy dependence of B. In order to take this into account we 
use the following form for the energy-dependent slope parameter 



B{W) = B(W ) +2a'jp \n(W/W f 



(3.6) 



with B(W ) = 4.0 GeV" 2 at an input scale of W = 40 GeV and following ZEUS we 
take a'jp = 0.1 GeV -2 . This clearly implies a reduction in the overall normalization 
of the cross section at large W, relative to the case in which a constant B is employed 
(for W = 200 GeV the reduction is about 14 % ). 

A comparison of the short-dashed and solid lines of fig. (H) show the significant 
effect of including the IV-dependent slope according to eq. (|3.6| ). The dotted curve 
shows the combined effect of running quark mass and B(W) on the solid curve (cf. 
also eq.( |3.5|) ). Both effects reduce the steepness of the energy dependence bringing 
the theory curves closer to the data. 

3.4 Skewedness 

Strictly speaking for exclusive processes we need to replace the ordinary gluon in 
eq. (|l.l| ) with the skewed gluon G(x\, 5, Q 2 ) (see p9| , and references therein, for a 
review of skewed parton distributions). For J/ip this should have a rather small 
effect, for T, with its much larger mass, it is certainly a large effect [[J. It is not 
obvious precisely how to sample the skewed gluon. Assuming collinearity of the 



10 



r-O 



0.03 



0.025 - 



0.015 




"1 1 

Fixed mass 
Running mass 



0.4 0.6 

b (fm) 



0.8 



Figure 3: The integrand multiplying a for fixed and running charm quark mass. The 
running mass effect is implemented for b < 0.3 fm and causes a suppression of the small b 
region. 



incoming and returning gluons (carrying momentum fractions X\ and — x^ of the 
proton, respectively) we have (x 1 P + q) 2 = M 2 5 , ((x 1 — x 2 ) P + q) 2 = M? 
photoproduction limit for the skewedness parameter, 5, this gives 
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(3.8) 

At first sight it appears there may be a danger of entering the ERBL |K| region 
(xi < 5, X2 < 0) in certain points in the phase space (in particular for symmetric, 
z ~ 1/2, configurations with small k 2 ). However, our ansatz protects us from this, in 
the case of J/ip, since we only use eq.( |l.l|) for b < 6qo, which corresponds roughly to 
k 2 > Q 2 Q . Even for an input scale of 1.0 GeV 2 this implies M 2 5 > 4 (m 2 + Q 



g) > 13.0 



GeV 2 . This is much bigger than the square of the J/ip mass: Ml(lS) = 9.59 GeV 2 . 



The 2S-state, ip', with a mass of M 2 ,(2S') 
which merits further investigation. 

For J/ip photoproduction we will use 



13.59 GeV 2 is a more marginal case, 
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11 



(cf. eq.( p7lf) with x' min = 5 = M v /W 2 ). This choice obviously guarantees that x' > 5 
for all b, so restricts us just to the DGLAP region (with this assumption we can never 
enter the ERBL region). With A = 10, m c = 1.5 GeV, (b) 2 = 10* (hc) 2 /4m 2 c » (0.21 
fm) 2 and this reduces to 

*<« = '< 1 + <i&- (3 ' 10) 

In our computer codes, the divergence at b — > 0, in the numerically unimportant 
very small 6-region, is regulated by hand by adding a very small number to b 2 in the 
denominator. The skewed gluon density, G(x',5,Q 2 ), is sampled at four-momentum 
scale Q 2 = X/b 2 ~ 0.39/(6 (fm)) 2 GeV 2 , hence some regulation of this scale is also 
implemented at very small b 2 (in practice we don't let the scale get larger than 
Q 2 = 100 GeV 2 ). Figs.f^J^) illustrate the 6-dependence of x'/S and Q 2 for different 
values of parameter A. 
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Figure 4: The effective momentum fraction at which the gluon is sampled, divided by 5, 
for several values of the scaling parameter A. 



In order to implement skewedness we use the evolution package developed by 
Freund and Guzey |HJ], which is based on the CTEQ code p2| . The skewed parton 
distributions (SPDs) are implemented in a systematic way underneath the integral 
in b. For a particular value of W, 5 is fixed. For this fixed 5, the skewed gluon is 
sampled on a grid, G(x', 5, Q 2 ') and the values saved to an array. The integration 
over b is performed numerically folding in the M^ 2 -independent piece from the wave- 
functions with a calculated by interpolating the skewed gluon array appropriately 
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Figure 5: The effective four- momentum scale at which the gluon is sampled for several 
values of the scaling parameter A. 



using a splines-based interpolation routine. Fig.(|GJ) shows the skewed gluon versus 
the conventional one at a scale typical for J /if; photoproduction for x' close to, but 
larger than, 8 = 9.6 x 10~ 4 (the value which corresponds to W = 100 GeV). Also 
shown is the SPD with very small S = 10 -7 which coincides with the conventional 
PDF, illustrating the limit G — > xg as 5 — > 0. 



Following closely the prescription for a given in [|lj] we impose the following 
behaviour at b = b n : 



a(b n ,W 2 ) = a(k 




(3.11) 



with a(b n , Wo) = 24 mb, e = 0.08 and Wq = 31 GeV. The latter value is chosen to 
coincide with the choice xq = 0.01 made in JT2|. A simple linear ansatz is used for 
the skewed a in the region b cr u < b < b n (we connect b cr u and b n with a straight 
line, b cr it is the point in b at which cr(b, W 2 ) = crib^, W 2 )/2 ). For small W 2 taming 
isn't required in the perturbative region, b < 6q , so the straight line starts at 6q - 
The dashed lines of fig. (0) show the effect of including skewedness in the dipole cross 
section relative to using the standard PDF (solid lines) at W = 100, 300 GeV. 

Close to the boundary at b = 6q = 0.39 fm, x' ~ 5 and Q 2 ~ Qq. If one 
assumes the standard input at the boundary for skewed evolution as we do, the 
difference between the use of skewed and standard distributions is minimal close 
to bq Q . For very small b, although there is a large evolution scale, Q 2 ^> Q%, the 
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Figure 6: Skewed and conventional (PDF) gluon distributions at a scale typical of J/ip 
photoproduction, Q 2 = 6.76 GeV 2 . The value of 5 = 0.00096 corresponds to W = 100 GeV. 
This figure explicitly illustrates that the skewed distribution reduces to the standard one 
in the limit of zero skewedness. 



gluon is sampled at x' ^> 5, so the overall effect of skewedness is expected to be 
fairly small everywhere. This is illustrated explicitly in fig.® which shows the 
integrand of eq.( [2.7|) using skewed and standard gluons at two different energies 
(W = 100, 300 GeV). The maximum effect, of about 10%, is seen close to the peak. 
For the skewed case it is interesting to show how the 6-distribution changes with 
energy, since this reflects the interplay of short and long distance contributions at 
different energies. One observes in fig.(|9|) that the peak shifts to the left as the energy 
increases and becomes more narrow, indicating an increase in the relative importance 
of short distance effects in this region. Examining such plots at very high energies 
reveals how the unitarity corrections begin to set in. One can start to see this in 
the shape of the curves at W = 300, 500 GeV, where the taming restriction begins 
to remove part of the distribution to the right of the peak. Although the taming 
corrections start to bite around 300 GeV they take a long time to tame the majority 
of the (fairly broad) peak in b. In practice this implies that we expect the energy 
dependence to be tamed very gradually to the soft Pomeron one. 



Finally, in fig. (10) we show the effect of including skewedness at the cross section 
level. The solid line corresponds to standard gluon PDF in a (with running mass 
and shrinkage implemented). The long-dashed line corresponds to replacing the 
ordinary gluon distribution with the skewed one. The overall effect, an enhancement 
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Figure 7: The dipole cross section a including skewedness (dashed curves) and without 
(solid curves) for two photon-proton energies. 



of approximately 10%, is seen strongest at high energies where the small dipoles play 
an increasingly important role (see fig.@). The shape in W still appears to be too 
steep. It is possible that if we allow unitarity corrections to play a role earlier we 
may be able to fit the high energy data better. 

3.5 Including the real part of the amplitude 

Having implemented the skewed gluon in the imaginary part of the amplitude we 
now reconstruct the real part using the analytic properties of the amplitude [|33| 



Numerically we achieve this by performing a two-power fit to Xm A(W) over a very 
wide range in W ( 10 < W < 900 GeV) using the form 



W 2 W 2 



the real part is then given by 



W 2 



Tie A(W) = ai (— f ) pi tan( 




WrV "" 2 



+ a 2 (ttto) tan (^-) • 



WrV "" 2 



The ratio of real to imaginary parts, (3 is then obviously given by 
tan(7rpi/2) + a 2 /a 1 (W 2 /H/ 2 ) P2 ~ P1 tan(vrp 2 /2) 



P 



1 + a 2 / ai {W 2 /Wi)P^-P^ 



(3.12) 



(3.13) 



(3.14) 
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Figure 8: A comparison of the 6-integrand of eq.( p.7| ) with (SPD) and without (PDF) 
skewedness at two different photon-proton energies. 



This increases the normalization of the cross section by j3 2 % (cf. eq.( |2.9| )). 



We achieve an excellent two-power fit using MINUIT |34| and get the following values 
for the fit parameters: 



ai = 0.00104, a 2 = 
Pl = 0.226, p 2 



-0.000794 
-0.0775 . 



(3.15) 



Fig. flTTD illustrates that j3 decreases as a function of W as the smaller power 
becomes less significant (see curves labelled A = 10). The short dashed line of fig.(|lO|) 
shows the overall energy-dependent enhancement (about 20% at high energies) when 
this implemented at the cross section level. This improves the description of the low 
energy data and in terms of shape, but leads to a worse overshoot at high energies. 

4. Decreasing A and predictions for THERA 

At large Q 2 , within the leading and next-to- leading logarithmic approximations of 
perturbative QCD, there exists a rather straightforward relationship between trans- 
verse size and the relevant four-momentum scale for the process concerned, i.e. that 
they are inversely proportional j|, [J. The constant of proportionality, A, was de- 
termined in [|| by an iterative averaging procedure involving the integral in b for 
the structure function F L : X — (bp L ) Q 2 - A value of A ~ 9 was obtained from this 
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method and this was found not to vary too much with x and Q 2 at small x ~ 10~ 3 , 
provided Q 2 > 10 GeV 2 . At larger b 2 , where applicability of perturbative QCD can 
not be justified, this relationship may break down, i.e. A may be found to depend 
on b 2 or equivalently on Q 2 . Indeed, in we observed that the value coming from 
this procedure was seen to deviate when a wider kinematic range was considered. 
During the work carried out for |12[] it was found that although A was changing, if 
kept within reasonable limits (4 < A < 15) the result for Fl only changed by a few 
percent (see fig. (5) of [I2| for an explanation of this approximate scaling). Having 
made this observation, a value of A = 10 was chosen for convenience (and because it 
most closely corresponded with average value of b in Fl) in the analysis of inclusive 
Structure Functions. 

At smaller Q 2 , we need to account for non-perturbative QCD effects, where the 
relationship between b 2 and Q 2 is not so straightforward (in particular, b 2 does not 
tend to infinity as Q 2 — > 0). To account for this slower dependence on Q 2 , for smaller 
Q 2 , and to test the sensitivity to non-perturbative QCD effects we diminish A at all 
b. It is logical that decreasing A will have a much bigger effect on processes such 
as J/ip photoproduction which contain greater contamination from large distances 
and are sensitive to scales at which the gluon density is changing rapidly at small x 
{Ql < Q 2 < 10 GeV 2 ). 
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Figure 10: The effect of including both skewedness and (5 at the cross section level. 
The solid curve uses the ordinary gluon PDF with running mass and W-dependent slope 
included. The dashed curves include skewedness. The short-dashed curve also includes the 
real part of the amplitude. 



Changing the value of A has several effects. Firstly, the position of the input 
scale, 6q , i n b shifts which affects a directly (recall that b n which also specifies a is 
fixed). It also directly influences the scales at which the gluon distribution is sampled, 
Q 2 and x', (cf. eqs.(|3.2|, |3T9| ) and figs.(f|,|5])). Since the light-cone wavefunctions do 
not depend on A (except implicitly though the Q 2 -dependence of m c>r ), this change 
has the effect of squeezing or dilating the perturbative region in b. Decreasing A 
decreases 6qo = vA/Qo and so diminishes the perturbative region almost without 
modifying a in the perturbative regime. The effect on a is explicitly illustrated in 
fig. ([13]) which shows the dipole cross section for several values of x and A = 4, 10 (for 
the purposes of this figure we neglect the 10% skewedness effect). This should make 
the cross section rise less steeply with energy as it enhances the non-perturbative 
piece. To test this we reran our code using the lowest value which left the results for 
F L unaffected, i.e. A = 4 \ 

The effect on the cross section was rather dramatic and is shown in fig.(|i"2|) (the 
real part of the amplitude, calculated as described above, is included, cf. fig.(|Tl"|)). 



^Incidentally and accidentally, this is the value advocated in the papers of Gotsman, Levin and 
Maor and collaborators for all Q 2 . Sec |35f for the latest version of their model which includes 
shadowing corrections. It is applied to J/ip production in pq|. 
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Figure 11: A plot to illustrate the relative size of the real part of the amplitude as 
a function of energy (for CTEQ4L partons). In the case A = 4 a two-power fit to 
MINUIT fit produces two positive terms with positive powers, hence (3 calculated from 
the fit increases slightly as a function of energy (the fit parameters in this case are: 
a ± = 0.000386, pi = 0.291, a 2 = 0.000186, p 2 = 0.0272). 



The cross section is increased in the fixed target region and suppressed in the high 
energy region. Both effects move it in the direction of the data. From ng.(|nj) we 
can see that big differences between A = 4 and A = 10 only emerge in the region 
around b = 0.3 fm. Hence, for 200 < W < 300 GeV, where the main contribution 
to the cross section originates from this region in b, the largest relative difference 
between the final cross sections is observed (see fig.(|l^)). The difference does not 
increase further with energy since the contribution from smaller b gradually becomes 
dominant (see fig.(|9|)). 

As an additional cross check we re-examined our description of F 2 (x, Q 2 ), using 
A = 4. The results for selected values of Q 2 are shown by the dashed lines in 
fig. (14). A comparison with the solid lines (A = 10) illustrates that F 2 is also fairly 
insensitive to this change in A and we still provide a reasonable description of the 
HERA data. Also shown (dotted curves) is the evaluation of F 2 using the latest 
MRST leading-order gluon distribution [|I7] , with A = 4. This appears to get closer 
to the data than the earlier CTEQ4L distributions (at least in the deep inelastic 
region). Taken together both the MRST and CTEQ4L curves illustrate the spread 
of predictions available from modern leading-log fits and hence give some indication 
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Figure 12: Decreasing the value of A from 10 to 4 improves the agreement with data 
dramatically for CTEQ4L partons. The short-dashed implements MRST input partons, 
with A = 4. 



of the theoretical uncertainty of the description. 

Extrapolating to higher energies we observe a rather broad spread of predictions 
for J/^-photoproduction in the THERA range in Fig. flT5|) . Since, for CTEQ4L 
partons, the case A = 4 does a better job on the lower energy data, we favour 
the dashed curve as our prediction for THERA. Despite the effect of the unitarity 
correction in the integrands, evident in fig.(|9]), the overall taming effect on the energy 
dependence at the cross section level appears to be rather mild within the considered 
model. This is because the contribution of very small b for which taming is still not 
important becomes more and more significant as the energy increases. In order to 
illustrate the sensitivity of the predictions to the choice of input parton density set 
we also show in Figs.( |T2"| , |T5| ) (short-dashed) curves for the latest MRST leading-order 
partons [|1(J. The latter assume an analytic form for the small-x behaviour of the 
input gluon which decreases as a function of l/x, which explains the milder energy 
dependence. 

Small A (e.g 4-5) appear from the J/^f analysis to be required where large b is 
significant. This is strongly related to the relative influence of perturbative and non- 
perturbative contributions. It seems that the current uncertainty associated with the 
small- x gluon distribution prevent any strong statement being made about A at small 
b. One may choose A in the range 5-15 and the results for physical quantities don't 
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Figure 13: The interaction cross-section, a for CTEQ4L, linear ansatz, for x = 
0.01,0.001,0.0001, A = 4,10. A value of Q 2 = Q\ = 2.56 GeV 2 is used in the ansatz 
for 6-dependent scales in a. 



change in quality (at present this applies to F 2 and hard vector meson production, 
since Fl is not yet measured). Making this change is perfectly allowed within the 
leading-log accuracy of the model (it corresponds to an attempt to mimic some NLO 
corrections). This observation shows that an estimate of the region where unitarity 
corrections become significant [EL [T2], |5P| will contain large uncertainties before a 
NLO calculation of cross section is made. 



5. Evaluation of a'p 

Since a'p is rather sensitive to the physics relevant for the taming of parton dis- 
tributions, in this section we develop a more sophisticated approach to calculating 
it. 

We would like to point out that the measurement of a'p recently reported by 
ZEUS [215] is entirely consistent with the results of our model. For large transverse 
sizes one expects a contribution of ajp(soft) = 0.25 GeV -2 , whereas for very small 
sizes a negligible contribution is expected. We can see this by introducing a simple, 
but reasonable, model based on dimension analysis: 



a'jpib) = 0.5 



6 2 + 6 2 



) GeV- 



(5.1) 
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Figure 14: A comparison of the inclusive cross section Fi using the dipole cross section 



with A = 4 (dashed line) and A = 10 (solid line), with a selection of the ZEUS [37| and HI 
1 38] data. Also shown (dotted curves) are the values obtained using leading-order MRST 
partons with A = 4. 



This model is designed to give a'pib = b n ) = cvj>(soft) = 0.25 GeV -2 , and to tend 
to zero quadratically at small b. There is a gradual decrease of a' F with increase of 
energy, as small b configurations become more important. We quantify this in fig.(|l^) 
by plotting the energy dependence of the average of a' P denned as (cf. eq. (|2.7| )) 



<OLjp> 



Jbdba'jp(b)I z (b)a 



(5.2) 



fbdbi z (b)a ' 

for A = 4, 10. We see that this model is in broad agreement with the ZEUS value 
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Figure 15: The photoproduction cross section for J/ifj in the THERA region using two 
different values for the scaling parameter A and CTEQ4L partons. The short-dashed curve 
shows the prediction using MRST leading order partons and A = 4. 



T-2 



13,11: a'jp = 0.098 ±0.035 ±0.05 GeV" 

Strictly speaking, the energy dependence apparent in fig.(|i"6"D may be slightly too 
strong when considering its effect on the t-slope parameter B. The average b, at a 
particular energy, is built in to the value of B(Wq) at the normalization point (cf. 
eq.(|3.6|)), but is expected to decrease with energy. We illustrate the point using the 
following phenomenological parameterisation for the slope of ^-dependence: 

(<b 2 >/4 + r 2 N ) 



B(W\Q 2 



+ 2a' ]P \nW 2 /W 2 



(5.3) 



Here <b 2 > is the average distance between c and c, denned in an analogous way to 
the average a'jp (cf. eq.( |5.2| )). In this phenomenological model a possible presence of 
ln& 2 /&Q terms is ignored. This formulae may be thought of as defining a convention 
for a new parameter, a', related directly to the logarithmic energy dependence of the 
slope parameter, B: 

OB 



a. 



a 



d (4 ln(W/W )) ' 
<b 2 (W)>- <b 2 {W )> 



+ a 



ip ■ 



(5.4) 



A8\nW/W 

However, it turns out that in HERA kinematics the difference between a' and a'jp 
from this effect is only about —0.01 GeV -2 . 
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Thus, the observed energy dependence of the average value of b in hard ampli- 
tudes (which are influenced by DGLAP evolution) leads to a more complicated and 
process-dependent energy dependence of the slope parameter, B, than that which 
is predicted by universal soft Pomeron exchange. At THERA energies and beyond, 
where the taming region starts to dominate, one expects that a' F should start to 
increase again. 
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Figure 16: The average shrinkage parameter, a' (using the simple model of eq.( |5.1| )) as a 
function of W (GeV), for two values of A. 



In our picture we expect both B and a'p to change with energy and photon 
virtuality as the balance between short and long distance contributions shifts. A 
dedicated forward detector for measuring scattered protons, such as the one recently 



proposed by the HI collaboration p0[ , would allow this issue to be investigated in 
detail. 



6. Discussion and open questions 

Using A = 4 and CTEQ4L partons our model still seems to overshoot the available 
data somewhat (cf. long-dashed line in fig.(|i"2|)), whereas MRST undershoot ((cf. 
short-dashed line). We would like to reiterate the point that the gluon distribution 
at small x ~ 10 -4 is not very well constrained by the current data (which mainly 
provides an indirect constraint via DGLAP-driven scaling violations of F2). This fact 
is reflected in the wide spread of numerical values for the gluon distribution at small 
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x in the currently available partons distributions " , of which we have shown only two 
here. Hence we are not too concerned by the fact that our curve appears to overshoot 



particularly the HI data [23| at the highest HERA energies using CTEQ4L. The 
framework we have described is general and can clearly use any leading-log parton 
density set that is available. However, we are encouraged by the improvement in 
agreement in the overall shape, which allows us to get closer to the data over a wide 
range in energy. 

It is natural to ask if higher order Fock states (\ccg> etc), which are only for- 
mally suppressed by a s , are important in the photoproduction of J ftp. The rationale 
employed in J3J was that using the solution to the Schrodinger equation (boosted to 
light cone) should take into account the most important corrections for large trans- 
verse distances. Recent studies of the effect of radiative corrections on the leptonic 



decay width |22| indicate that perturbative corrections for small distances may well 
be large. At present it is unclear what effect this will have on the photoproduction 
of J/ip. From the theoretical side what is required is a complete next-to-leading- 
log calculation. This question could be addressed phenomenologically by examining 
the relative characteristics of various diffractive charm measurements (exclusive J/ip, 
J/i/j + 1 jet, and open charm). The question is even more urgent in the case of i/j' 
which is obviously a larger bound state and contains a node. 



7. Conclusions 

We have investigated the photoproduction of J/ip in the context of the QCD-improved 
dipole model introduced in fll2"| . This model directly incorporates a contribution 
from long distances which is responsible for the low energy production. As the en- 
ergy increases, the short distance mechanism of perturbative two-gluon exchange 
becomes increasingly important. Overall our description of the data over the whole 
measured range is improved relative to analyses which only take the perturbative 
contribution into account. At very high energies we incorporate taming or uni- 
tarity corrections and present predictions for the cross section in the THERA range 
(250 < W < 900 GeV). In our model, unitarity corrections affect smaller and smaller 
transverse sizes at progressively higher energies. For J/ip photoproduction, which is 
sensitive to a broad range in transverse sizes, it turns out that in the THERA region 
there is always a significant contribution to the cross section that rise quickly with 
energy. This implies that the taming of the growth of the cross section with energy 
is rather gradual (the precise details depend on the choice of input density). 

We show that the skewedness of the amplitude induces a relatively small en- 
hancement in the cross section (approximately 10%). This effect is likely to be 



See http://durpdg.dur.ac.uk/HEPDATA/PDF3.htm] 
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swamped by other uncertainties associated with the ^-dependence and the light-cone 
wavefunction of the vector meson. 

Of much greater numerical importance in describing the available data is to 
modify the balance between short and long distances contributions. In our model, 
this may be controlled by changing the scaling parameter A which sets the relationship 
between four-momentum scales and transverse sizes in the interaction cross section. 
For CTEQ4L partons, the shape in energy of the J /if) photoproduction data appears 
to favour a change to a lower value (A = 4) for large b than that which was derived in 
[|], H, |12| for small b (A = 10). The quality of the description of the structure function 
F2 is relatively unchanged by this modification. The largest uncertainty comes from 
poor knowledge of the numerical size of the leading-log small-x gluon distribution. 
This leads to a rather broad band of predictions for J /if) photoproduction in THERA 
region. Studies of the photo- and electroproduction of other vector mesons, within 
the same framework, is underway. 
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